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Abstract 

We propose a new mechanism to block the active-sterile neutrino oscillations 
in the Early Universe. We show that a typical consequence of theories where 
the lepton number is spontaneously broken is the existence of a coherent 
cosmological Majoron field with a strength proportional to the lepton and 
baryon numbers of the Universe. This field interacts with leptons and changes 
the potentials relevant for neutrino oscillations. If the scale of lepton number 
symmetry breaking is of the order of 1 GeV then a Majoron field and lepton 
number asymmetry of the order of the baryon asymmetry are strong enough 
to block the active-sterile neutrino oscillations with the atmospheric neutrino 
mass gap which otherwise would bring the sterile neutrino into equilibrium 
at the big bang nucleosynthesis epoch. 

PACS numbers: 14.80.Mz, 14.60.Pq, 95.30.Cq, 98.80.Ft 
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I. INTRODUCTION 



The explanation of the present neutrino puzzles M may require the existence of one or 
several species of extra light sterile neutrinos 0. In particular, the sterile neutrino could be 
relevant for the explanation of the atmospheric neutrino problem (ANP) || in the presence 
of a significant — v s mixingQ. The typical required values are 5m 2 &tra ~ 3 x 1CT 3 eV 2 and 
large mixing angle, sin 2 2# atm ~ 1. 

On the other hand, for such a parameter range one can encounter a contradiction with 
the Big Bang nucleosynthesis (BBN) bounds |J on the number of extra light particle species; 
namely, according to the analyses of ref. J?J, the sterile neutrino comes into equilibrium with 
the particle thermal bath via — v s oscillation, unless the condition 5m 2 sin 4 29 < 3 x 10 -6 
eV 2 is satisfied (updated constraints are given in ref. M for small mass differences, 5m 2 < 
10~ 7 eV 2 ), which is certainly out of the range of parameters needed to explain the ANP. 

However, it was found || that the — v s oscillations are suppressed at temperatures 
T < 3 MeV (the decoupling temperature of ^ r ) if the lepton number asymmetry at these 
temperatures is very high, namely L a > 10~ 5 (lepton number to photon number ratio). But 
this is 4 — 5 orders of magnitude larger than the observed baryon asymmetry of the universe 
(B < 10 -9 ) and in the most generic baryogenesis context one can expect that L ~ B (e.g., 
in the context of grand unified theory (GUT) baryogenesis or leptogenesis PH|JTT[] this is 
because the B + L non- conserving sphaleron processes redistribute B and L among each 
other P2]-4T1[). The same is true in the context of the electroweak baryogenesis (B — L = 0). 



In the Affleck-Dine mechanism B and L can in principle be independent of each other, but 
still of the same order. 

It has been shown that at much lower temperatures (T < 100 GeV) neutrino 

oscillations can actually produce a rapid increase of the lepton asymmetries from the initial 
very small values up to the order of 0.1. This however only occurs for negative 5m 2 cos 26 
and very small active-sterile mixing angles, not directly relevant for the ANP. 

In this paper we show that a lepton asymmetry as small as the present baryon asymmetry 
may be enough to block the sterile neutrino oscillations. The necessary new ingredient is the 
existence of a coherent Majoron field in the Early Universe and a low scale of spontaneous 
breaking of lepton number, Fl ~ 1 GeV (The Majoron [|Hj] is the massless Nambu-Goldstone 
boson in models where the total or any partial lepton number is spontaneously broken.). 

While it has been common wisdom that due to their derivative coupling nature [17] 
Nambu-Goldstone bosons cannot mediate long range interactions, it has recently been 
demonstrated [|18|] that a coherent source of a Majoron field, to be specific, is formed when- 
ever the corresponding broken lepton number suffers a net increase or decrease in a certain 
region of space. The processes that violate this lepton number can be the very neutrino os- 
cillations as exemplified in previous papers |18| or any other reactions. In the present work 
we show that a Majoron field can be produced due to lepto- and baryogenesis processes in 



1 The recent Super-Kamiokande data can be explained by — v T oscillations while the situation 
where the ANP is exclusively due to v^ — Vs oscillations is disfavored Q. However, the more general 
case where oscillates into v T and v s with comparable rates is completely consistent with the 
data f|. 
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the Early Universe. The Majoron field interacts with neutrinos with a strength inversely 
proportional to the lepton breaking scale, Fl. If Fl is around 1 GeV, a lepton asymmetry as 
small as L ~ B ~ 10~ 9 can block the oscillation into sterile neutrinos with Sm 2 ~ 3 x 10~ 3 
eV 2 no matter how large their mixing angle is. That is our thesis. 

The origin of the Majoron field is elaborated in section [V] and its role on neutrino oscilla- 
tions into a sterile neutrino in section |VI|. But first we build in section |T| a specific model of 
neutrino masses with spontaneous breaking of lepton number and in section [TV] we derive the 
relations between the particle asymmetries in the Early Universe and the present baryon 
number. The model aims to fit the present known observations from solar, atmospheric 
and terrestrial neutrino experiments, including the Liquid Scintillation Neutrino Detector 
(LSND) result [19|, with an extra sterile neutrino. However, it should be emphasized that 
the mechanism we propose of suppression of oscillations into a sterile neutrino at the BBN 
epoch, based on the existence of a Majoron field, does not depend on the particular model 
or set of neutrino mass parameters. The only fundamental assumptions are that the lep- 
ton number (or a partial lepton number) is spontaneously broken and the breaking scale 
is around the 1 GeV magnitude. In the next section we also make the point that in the 
absence of the LSND neutrino mass gap the oscillations of atmospheric and solar neutrinos 
into a sterile neutrino are no longer correlated with each other. In the last section we draw 
our conclusions. 



II. NEUTRINO MASSES AND MIXING 

The existence of a fourth, sterile, neutrino has been suggested as it is the only way of 
reconciling the atmospheric, solar, and LSND neutrino oscillation evidence and their very 
different Sm 2 mass gap scales. The wide mass gap (0(1 eV)) that is necessary to explain 
the LSND result in terms of v e — mixing requires that the neutrino mass pattern should 
have a two-doublet structure |20| : one of the doublets consists of v e and v s or z/ r , or a 



linear combination of both, and is responsible for the solar neutrino deficit, and the other 
one, responsible for the atmospheric neutrino anomaly, consists of v ^ and v T (or a linear 
combination of v s and u T ). These doublets are separated by the LSND mass gap. 

It has been shown pT] . |22| that even after the recent SNO observations [^| both the 
sterile and active neutrino oscillations are viable solutions of the solar neutrino problem as 
well as a more general superposition of both. On the other hand, the atmospheric neutrino 
data seem [|j] to favor the v T solution against the sterile neutrino case but the analysis 
of the most recent data still allows a quite large relative probability, more than 50%, of 
oscillation into v s (the larger the probability, sin 2 £, the smaller the allowed Sm 2 range). 

A consequence of the LSND large mass gap and the limits from reactor disappearance 



experiments such as Chooz |24[ is that the solar neutrinos v e and the atmospheric neutrinos 
must oscillate into states that are essentially orthogonal to each other. In other words, if 
the solar electron neutrinos oscillate into the linear combination v s = cos£ v s — sin£ u T , then, 
the atmospheric muon neutrinos necessarily oscillate into the state Up = sin£ v s + cos£ v T . 
However, the situation would be totally different if the LSND evidence was not present. 

To be more specific, let V\, V2 and 6 & be the mass eigenstates and mixing angle responsible 
for the solar neutrino deficit and z/ 3 , i/ 4 and 8 atm the states and mixing angle relevant for 
atmospheric neutrinos. The two pairs are separated by the LSND mass gap and no other 
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specific mass hierarchy has to be assumed. The reactor experiments constrain the mixing 
matrix elements U e3 , U eil U^i, and to be small [pOfl , but not U Ti or U si . If one neglects 
all the mixing angles that are necessarily small and irrelevant to explain the present bulk of 



data, the mixing matrix is given as 

V\ = cos9 Q v e — sin# Q (cos£ v s — sin£ z/ T ) , (la) 

z/ 2 = sin 6*0 v e + cos6 , (cos^ v s — sin£ v T ) , (lb) 

z/ 3 = cos6> atm v ^ - sin # atm (sin £ v s + cos£ v T ) , (lc) 

1/4 = sin # atm z/ M + cos # atm (sin £ i/ s + cos f z/ T ) . (Id) 



The mass eigenstates z/ x and u 2 are separated by the gap 5m 2 & and z/ 3 , z/ 4 by Sm 2 ^. 

Clearly, the less the atmospheric neutrinos oscillate into the sterile neutrino the more the 
solar neutrinos have to oscillate into u s . There is a potential clash in the future if both solar 
and atmospheric neutrino experiments happen to constrain the respective sterile neutrino 
solutions to less than a 50% probability. In that case the conflict with the LSND data will 
be insoluble, which will call for new results from MiniBooNE | ]25[| , the next new independent 
accelerator experiment. Suppose for a moment that the LSND evidence does not exist or is 
going to be ruled out by the MiniBooNE experiment. We would like to stress that this does 
not rule out the sterile neutrino as a possible protagonist in the other solar and atmospheric 
neutrino problems, not even if both of them exclude dominant sterile neutrino solutions. 
On the contrary, the absence of the LSND mass gap increases the freedom in the neutrino 
mixing parameters. Then, the two-doublet mass pattern is no longer inevitable and the role 
of v s in the solar neutrino deficit is completely decoupled from its role in the atmospheric 
neutrino oscillations. 

As a matter of proof we make explicit an extreme case, namely, where the atmospheric 
neutrinos oscillate into v T or v s with arbitrary relative probabilities, while the solar neutrinos 
oscillate exclusively to v T and but not to u s . The mixing matrix can be described as 
follows: 

ui = cos 9 Q v e — sin^ Q (— sin a + cos a z/ T ) , (2a) 
v 2 = sin 6*0 v e + cos# (— sin a + cos a u T ) , (2b) 
z/ 3 = sin (3 v s + cos (5 (cos + sin a z/ T ) , (2c) 
z/ 4 = cos P v s — sin (3 (cos a + sin a v T ) . (2d) 

As far as the mass spectrum is concerned, the mass eigenstate u 3 is separated from the other 
three by mass gaps that are in the atmospheric neutrino range ~ <5m^ tm ~ 3 x 10~ 3 eV 2 . 
v\ and z/2 are almost degenerate and separated by the solar neutrino mass gap and finally, 

is only subject to the condition m\ — m\ ~ Sm 2 ^^ as it, like v 3 , does not participate in 
the solar neutrino oscillations. The mixing angles relate to the atmospheric mixing angle 
as cos # atm = cos a cos f3. The atmospheric neutrinos oscillate into v T with a probability 
proportional to sin 2 a cos 2 (3 whereas the probability of oscillation into v s is proportional to 
sin (3. The ratio between them is given by tan 2 £ = tan 2 /?/ sin 2 a. It is clear that the solar 
neutrinos do not oscillate into v s . 

This just shows that the potential problem raised by the possibility that the atmospheric 
neutrinos oscillate significantly into a sterile neutrino with its consequences for BBN is not 
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necessarily linked to the solar neutrino solutions and does not depend on the LSND evidence 
although it has been motivated by the coexistence of all three kinds of observation. In the 
present work we want to present a solution and a mechanism to block the oscillations of 
muon neutrinos into sterile neutrinos in the Early Universe at the time of BBN. The idea 
does not depend crucially on the particular neutrino mixing pattern but the actual numbers 
vary, of course, from model to model. We worked out in detail a particular model that is 
suitable to encompass all three types of neutrino oscillation evidence, including LSND. 



III. NEUTRINO MASS MODEL 



The seesaw mechanism [26] can be incorporated within a model where the lepton number 
is spontaneously broken at a relatively low energy scale by adding to the standard lepton 
doublets £; and charged singlets two heavy sterile neutrinos per lepton generation, N Li 
and Ng (left-handed), with lepton numbers +1 and —1, respectively. The additional light 
sterile neutrino, v s (left-handed), has lepton number L s = —3. The most general Yukawa 
interaction Lagrangian in the lepton sector is written in Majorana matrix form as 



c Y = i^ T aw + H.c. 



(3) 



where ^ = (ef , £ h u s , N Li , Ng 



and A4 is the symmetric matrix 
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(4) 



The omitted elements are obtained by symmetrization. Hi and H 2 are two standard Higgs 
doublets under SU(2) and a is the singlet scalar field with lepton number L a = — 2. h s is a 
3x1 column and h e , h^, hi, Hr, and M are 3x3 matrices. 

Before lepton number spontaneous breaking the heavy sterile neutrinos form Dirac par- 
ticles, namely, Ni = Nn + Nm, with lepton number equal to 1 and masses Mi in the basis 
where M is diagonal: M = diag(Mj). After lepton and gauge symmetry breaking the light 
neutrinos acquire masses and mix with the sterile neutrino in a 4 x 4 Majorana mass matrix. 
Denoting the 3x3 active, 3x1 active-sterile and lxl sterile neutrino blocks respectively, 
as m vv , m us and m ss , we obtain in leading order in any basis where M is a real matrix 



m vv = /i JV M~ 1 /i L (/i 7V M" 1 ) T (a) v. 
m vs = -h N M~ l h s (a) v 2 
m ss 



2 

2 • 



(M^hsfhaM-X (a)' 



(5a) 
(5b) 
(5c) 



where v 2 = (H 2 ). We take as reference scales m uu ~ 0.05 eV to account for the atmospheric 
neutrino anomaly and m us ~ 1 eV for the LSND v^ipy) —>■ v e {v e ) evidence. Since we also 
assume (a) ~ 1 GeV the element m ss ~ 10 -13 eV is completely negligible. 
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IV. ASYMMETRIES IN THE EARLY UNIVERSE 



At temperatures below the heavy neutrino masses M; > 10 6 GeV, the Dirac masses Mj 
still mediate scattering processes capable of producing the light singlet particles v s and a 
like £H 2 — > v s a. They can be studied in terms of the effective operators 

C eS = iiHtt^ijHi a + £ t H 2 ^u s a* + u s ~P^-u a a* 3 + H.C. , (6) 
2(a)v 2 (a)v 2 2(a) 3 v 2 

which also give rise to the light neutrino masses after spontaneous breaking of lepton number. 
One obtains the cm. cross sections of the scattering processes 1) £i£j — > H 2 H 2 a, 2) aH 2 — > 
£iU s , and 3) era — > au s u s as 

(8vr) 3 {<t) 2 v% ((t) 2 ^ 

^ = ^r§4 « TvVt x ID" 2 eV 2 , (8) 
8% {a) 2 v 2 (a) 2 v 2 

a 3 = -^^^ » x 10- 28 eV 2 , (9) 
(8vr) 3 ((j) 6 (ct) 6 v 1 

respectively, where we have summed over initial and final weak isospin states (y/s is the cm. 
energy) . 

In each case one compares the rate of collisions per particle, T = an ss 0.1 a T 3 (the 
boson number density is ~ 0.122 T 3 and the fermion number density rif ~ 0.091 T 3 ) with 
the Hubble rate /J ~ T 2 /10 18 GeV, assuming a total number of degrees of freedom around 
100. The scalar singlet a is produced through the processes 1) £i£j — > H 2 H 2 a, £iH 2 — * £jH 2 a, 
and H 2 H 2 — > £j£j<r with cross sections o"i, 2<7i/3, and a"i/3, respectively (for i 7^ j), which 
gives a total rate per Hubble time r^tf" 1 w 2(w/w2) 4 T 3 /10 15 GeV 3 (u ~ 174 GeV is the 
electroweak breaking scale). This shows that a is in thermal equilibrium at temperatures 
larger than T a ss 10 5 GeV, if one takes v 2 = v. 

The sterile neutrino is produced in the processes 2) aH 2 — * £iu S) £i<y — > -f^s, and 
£j^ 2 - ► with cross sections o" 2 , <J 2 /2, and (J 2 /2, respectively, and a total rate Y S H~ X ss 
(v/v 2 ) 2 T/T s , which makes the decoupling temperature of the light sterile neutrino T s ss 
4 x 10 6 GeV. If one or more heavy neutrinos iVj have masses under that value, v s may 
decouple when some of the iVj degrees of freedom are still present in the Universe (note that 
Mi > 10 6 GeV). Finally, processes like go — > eri/ s i/ s are too weak to be relevant. 

Above T s the sterile neutrino and scalar singlet a are in chemical equilibrium with the 
lepton and Higgs doublets and their number asymmetries are constrained by the equations 
of detailed balance. The precise relations between the particle asymmetries depend on 
which particles and processes are in thermodynamical equilibrium at a given time. To 
be definite we assume that by the time the sterile neutrino decouples, B and L are only 
violated by electroweak instanton processes while B — L is conserved. On the other hand 
the right-handed electrons e# are not yet in chemical equilibrium and the quarks ur, 
may or may not be in equilibrium depending on the exact values of their Yukawa couplings 
and temperature T s . In either case the equations of detailed balance yield the particle 
asymmetries as functions of the B — L asymmetry. 
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At temperatures above T s the operators of eq. @ yield the chemical potential constraints 

fi a + 2m + 2fi H = , (10a) 
Us - Ha + H + (J>h = . (10b) 



The other constraints come from standard model reactions [T^JT^] • To be definite we assume 
that ur and cLr are in equilibrium at T s (the temperature at which they come into equilibrium 
increases with increasing Yukawa couplings and therefore with increasing number of Higgs 
doublets). The electroweak and QCD instantons and the Yukawa interactions imply that 

3/ig + fit = , (11a) 
2yU g - fi u ~ = , (lib) 
H Q - ~ He = , (11c) 

- Hu + = , (Hd) 

Ht - fi T -/iff = , (lie) 

where fi q , /ig, and /iff designate the flavor universal chemical potentials of the quark, lepton, 
and Higgs doublets, respectively, /t M , /^ those of the right-handed quark isosinglets and fi T 
the common chemical potential of the lepton isosinglets (Ir and tr. Since the electron singlet 
cr is not in chemical equilibrium its chemical potential /t e is an independent variable. We 
may assume that a baryon asymmetry originally produced in a GUT baryogenesis scenario 
is later communicated through electroweak instantons to the lepton sector (T < 10 12 GeV) 
but not to cr. In that case fi e remains zero until the cr Yukawa interactions come into 
equilibrium at temperatures lower than T s . 
A vanishing weak hypercharge implies 

3(/i 9 + 2/i u - fj, d - n £ ) - 2/i T - ji e + 2n H fi H = . (12) 

Here hr is the total number of Higgs doublets; hr = 1 if Hi and H 2 are the same field 
and ur = 2 otherwise. The above constraints and the condition /t e = leave only one 
independent variable. It is convenient to choose this as fi s because the v s abundance and 
number asymmetry are conserved after its decoupling. The other quantity that is conserved 
is B — L. Denoting the baryon number density as dB/dV = BT 3 /Q and likewise for L and 
B — L, one has 

B = Qfi q + 3fi u + fi d , (13a) 
L = - 3/t s , (13b) 
l>4 = 6fJ>i + 2/t T + /i e - 4/t CT + 2n N fx N , (13c) 

where L<f stands for the lepton number of all particles except v s and is the number of 
relativistic Dirac heavy neutrinos at a given moment. The decay processes iVj — > £jH 2 set 
the equation fi N = /x e + fi H . Putting everything together, one obtains 

B-L 4 =\U + 2n N - (Ns ~ Ns) ' (14) 

B-L (9 + nff)(17 + 2njv) - 15 - 4n H 



B-L 4 (9 + nff)(8 + 2njv) - 15-4nff 



(15) 
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After u a decoupling, B — L, B — Lj, and the v s number asymmetry N s — N s are all conserved. 
Although the above relations are strictly valid only when v s is in thermal equilibrium, one 
expects that the decoupling process does not introduce very large perturbations and one 
may use these results as a first approximation. They give the v s number asymmetry and 
B — as functions of the primordial B — L and number of heavy neutrinos that are 
relativistic when v s decouples. 

The next transition is the decoupling of the scalar singlet o at a temperature T a around 
10 5 GeV. This is close to the e# coupling epoch which starts at T e ~ (v/vi) 2 x 10 4 GeV. This 
temperature rises with increasing electron Yukawa coupling and if one assumes the existence 
of two Higgs doublets and V\ < v/3 (v ~ 174 GeV) then e# is already in equilibrium when 
a decouples. To be definite we assume so. One repeats the exercise with Eqs. (|Tl~D and 
(|I~2|), complemented with /i r = /i e and Eq. ( |10a| ), to obtain all chemical potentials in terms 
of \x a . Then Eqs. ( |13a|) and (|13c|) with = yield the relation between the a number 
asymmetry and the baryon and lepton numbers. Denoting the total lepton number of the 
standard model particles as Lg (Lg = L e + L M + L T ) one derives for nn = 2, 

12 (B-L t ), (16) 



N„- 




N s - 


N-s 


B- 
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B- 





23 

47 33 



23 65 + 22^ 
47 164 + 22n N 

23 65 + 22n 



(B - Li) , (17) 

(18) 



N 

where in the last two equations we used Eqs. ( |i~4"D and (|15D , keeping in mind that is the 
number of Dirac neutrinos iVj that are relativistic when v s decouples. Again, one expects 
that the above results remain a reasonable approximation when a decouples. From then on, 
B — Li, B — L, and the a and v s abundances are conserved by all effective interactions. 

When the temperature drops down to the electroweak phase transition the weak isospin 
and hypercharge are no longer conserved, contrary to the electric charge. The quarks and 
charged leptons form Dirac mass eigenstates whose well defined chemical potentials are 



subject to a new set of constraints |1J|, together with the neutrinos and charged Higgs and 
W bosons as follows: 

HW+ = Hu- Hd = Hv- He = VH+ , (19) 

3/i u + 3fi d + fi u + fi e = . (20) 

The latter constraint is due to sphaleron processes. On the other hand, the net electric 
charge is zero: 

3(4/i M - 2fi d - 2/i e ) + 2(2 + n H )nw+ = . (21) 

Definite B and L numbers are predicted in terms of the preexisting B — Li number, 
namely | |14| , 



B= 32 + 4UH (B — Li) , (22) 
98 + 13n H y ' ' K ' 

' e6 + 9 "»-<B-L t ). (23) 



98 + 13n H 



S 



These relations are preserved during the phase transition if there is no intrinsic electroweak 
baryogenesis. After that the sphaleron processes stop being effective and the baryon number 
is separately conserved. This allows us to predict the scalar a particle and sterile neutrino 
asymmetries in terms of the present baryon number. From Eqs. (|T6| ) and (|17|) , valid for 
two Higgs doublets (nn = 2), one derives the asymmetries and lepton numbers L( CT ) = 
—2(N a — N 5 ) and L( a ) = — 3(N S — N s ) carried by a and u s as 

24 31 

L M = -KB = ---B , (24) 

= ,£31 99 

(j 23 10 65 + 22njv 

It is important to notice that, after spontaneous breaking of lepton number, the B — L 
violating processes are too weak to be in equilibrium, in particular during the electroweak 
phase transition if it occurs after L breaking. If that was not the case the baryon and lepton 
numbers would be washed out. As soon as the sphalerons decouple the baryon and lepton 
numbers start to be separately conserved. This is also true after L spontaneous breaking 
because the L violating reactions are weak. The lepton number may only be significantly 
violated at much lower temperatures of the order of 1 to 10 MeV when neutrino oscillations 
from active to sterile neutrinos become possible. Another point is that, after L breaking, 
the lepton number Li a \ carried by the scalar singlet a still exists but is then associated with 
a coherent Majoron field. This is the subject of the next section. 



V. MAJORON FIELD 

The Majoron equation of motion is determined by the equation of conservation of the 
lepton number Noether current. The lepton current of the scalar field a with lepton number 
L a = -2 is 

J% = L a i (a* W - a W*) . (26) 

At the classical level the total lepton number of charged leptons, neutrinos and scalar a is 
conserved but electroweak instanton effects break L explicitly. B — L remains conserved and 
its equation of conservation reads as 

V^ + V M j; = 0, (27) 

where Jj is the L — B current of all the other particles, in our case leptons and quarks: 

J f = - - L f) ( n f - n i) v " > ( 28 ) 

where rif and nj are the particle and antiparticle densities and f M = (1, v) the macroscopic 
velocity vector. 

Before spontaneous breaking of the lepton number the a current is related to the a 
particle asymmetry, J^ 1 = L a (n a — n„) v^, but after lepton symmetry breaking, the mass 
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eigenstates are no longer the complex field a but rather the massive Higgs particle p and 
massless Majoron boson ip. They relate to each other as 

a = -^=(v a + p) exp(-z p/v a ) (29) 



and the lepton current is expressed as 



2 

P " 



JZ = L a v a ([l + ^) W>- (30) 



As emphasized in references |I8yi7| , after symmetry breaking the global symmetry is real- 
ized as an invariance under translations of the Majoron field whose equation of motion is 
determined by the still valid Eq. ( |2"7| ) of lepton number conservation. 

After lepton breaking the current J£ can be realized only through a coherent Majoron 
field (p; in other words, the expectation value of a has a variable phase: 

v 

W = ^expH^) . (31) 

We obtain for the current 

J£ = L uVu (l + ^Pj VV , (32) 

where the term (p 2 ) is an average over quantum fluctuations, i.e., the thermal bath of 
massive Higgs particles p. When the temperature of the p bosons, lower than the photon 
temperature when the number of relativistic degrees of freedom drops down one order of 
magnitude, is much smaller than the breaking scale v a , the (p 2 ) term can be neglected. 
Then, J£ = L a v a V^ip. In a homogeneous and isotropic Universe f depends only on time 
and 

J° = L a v a <p = F L <p . (33) 

The value of (p is subject to the equation of conservation (P7|). Integrating over space, the 
a lepton charge L( CT ) = J dV J® is determined at a given time by its initial value and the 
variation of the B — L number carried by leptons and quarks: 

L (a) (t) = L {a) {U) + (B- L) f (t) -{B- L) f {U) . (34) 

The initial value of L( CT ) is the lepton charge carried by the complex bosons a before spon- 
taneous lepton breaking. This value is proportional to the initial B — L or to the present 
baryon number, as Eq. ([24]) shows for the particular model we worked out. On the other 
hand, B — Lis possibly violated by neutrino oscillations only at very low temperatures when 
B is conserved. As a result, Lr a )(t) = —A a B — AL^ and the Majoron time derivative is 
obtained from Eq. ([33]) as 

if = --^(A a B + AL) (35) 
in terms of the baryon number and lepton number variation per photon 

B = — , AL = . (36) 

Notice that Lf and AL count only the fermion particles, charged leptons and neutrinos, but 
not the a field. 
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VI. NEUTRINO OSCILLATIONS 



Neutrino oscillations [g7| are governed by the neutrino masses and mixing angles and 
interactions with the background medium E3L which in turn depend on the temperature 
and particle number asymmetries. As far as standard model weak interactions are concerned 
the electron, proton, and neutron asymmetries, closely related to the baryon asymmetry, are 
too small to play a role in the oscillation of active neutrinos into sterile neutrinos. However, 
the neutrino asymmetries can in principle be much larger than the baryon asymmetry. 
Normalizing them to the photon density as 

K = , (37) 

the potential of the flavor v a = u e , u^, v T induced by electroweak interactions is given at low 
temperatures (T < <C M w ) by |2£| 

Vfew = ±V2G F n 7 (L a + L e + + L T t A a T 2 M^ 2 ) , (38) 

where A e = 55 and A^ T = 15.3 (the electron and nucleon asymmetries are neglected). The 
upper sign holds for neutrinos and the lower sign for antineutrinos. The sterile neutrino has 
no standard model potential by definition. 

In the case of interest, Sm 2 ~ 3 x 10~ 3 eV 2 , the thermal contribution proportional to 
n 7 T 2 prevents the oscillation into a sterile neutrino at temperatures above ~ 10 MeV. At 
smaller temperatures that term becomes ineffective and the active neutrino oscillates into 
the sterile flavor as in vacuum, violating the bounds on the number of light degrees of 
freedom at BBN [0]. Foot and Volkas || pointed out that if there is an initial asymmetry 
L a = L a + L e + L M + L T larger than 7 x 10~ 5 and 5m 2 < 10~ 2 eV 2 , the active neutrino 
v a cannot significantly oscillate into the sterile v s and the initial neutrino asymmetries are 
preserved until the active neutrino v a decouples or, in the case of u e , until the protons and 
neutrons stop being in equilibrium. In these conditions the BBN bounds on the extra light 
degrees of freedom are satisfied. These straightforward considerations have one price, which 
is the assumption of an initial neutrino asymmetry five orders of magnitude larger than the 
baryon asymmetry. f\ 

The situation changes if there is a Majoron field. Majorons, like any Nambu-Goldstone 
boson have only derivative couplings. As a result, a coherent Majoron field produces neutrino 
potentials proportional to its gradient []18[| . If A is the spontaneously broken lepton number, 
in general, any combination of partial lepton numbers, and A a the quantum number of the 
flavor z/ a , a Majoron field tp produces the potential 

V A = -—AaV^cp (39) 



2 In ref. Foot and Volkas explored the case where v T — > v s oscillations with —5m? > 10 eV 2 
and sin 2 26 < 10~ 5 create a lepton asymmetry large enough to block — > v s oscillations with 
ANP parameters. In any case, this cannot be a generic situation and requires some agreement in 
the parameter space - the masses and mixing of all neutrino species. 
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for the neutrino v a and the symmetric one for the antineutrino i> a , where = (l,v) is 
the neutrino four- velocity (|v| = 1 in leading order). In the present case the Majoron is 
associated with the total lepton number L and the Majoron field is a uniform field in the 
Early Universe given by Eq. (]35f) . Hence, it induces the potentials^ 

V L = F£ 2 n,L a (A a B + AL) , (40) 

where A a is a model dependent coefficient of order 1; A a ~ 3.2 in the case we are considering. 
The quantum numbers are L a = +1 (—1) for an active neutrino (antineutrino) and L a = —3 
(+3) for the sterile neutrino v s (z/ s ). 

The variation of the lepton number in the neutrino sector can be caused only by oscil- 
lations into the sterile neutrino because this is the only one with lepton number different 
from ±1 and the types of oscillations we are considering conserve chirality. The oscillations 
v a *^ v s [pa *^ v s ) produce a lepton number variation AL = AN a — 3AN S = —AAN S 
(AL = 4AiV 5 ); hence, 

AN S - AN S , . 

AL = -4 1 . (41) 

Combining Eqs. fl3"8|) and (|40|) , the difference between active and sterile neutrino potentials 
is 

V a -V s = V2G Fni (±L a - A a T 2 M^ 2 ) ± AF^n^A.B + AL) , (42) 

where L a = L a + L e + L M + L T and the lower signs apply to antineutrinos. It is now clear 
that if, in the case of standard weak interactions, an asymmetry L a > 7 x 10 -5 is enough 
to block the v a «-> u s , v a <-> u s oscillations then, in the presence of a Majoron field, the 
known baryon asymmetry B = (4 — 7) x 10 -10 can do the same if the scale of lepton number 
breaking Fl obeys the condition 

F^<^ixl0 5 = 5-9 GeV 2 . (43) 
7v 2Gp 

Such a low scale of lepton number breaking is perfectly consistent with the exist- 
ing bounds for this kind of singlet Majoron model [16,|32|,18| including the astrophysical 



bounds |I3|. The reason is that in scattering processes the Majorons couple primarily to 
neutrinos with strengths proportional to the neutrino masses g ~ i7i u /Fl, which are therefore 
negligibly small for the assumed neutrino mass spectrum even if Fi ~ 1 GeV. 



VII. CONCLUSIONS 

It is well known that an explanation of the atmospheric neutrino problem || in terms of 
oscillations of the muon neutrino into a, at least in part, sterile neutrino is in contradiction 



3 They should not be confused with the potentials induced by a thermal bath of Majoron particles, 
proportional to the neutrino masses, which could be significant for rather large neutrino masses 
like 17 keV as considered in ref. [31]. 
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with the BBN limits on the number of extra light degrees of freedom || . Indeed, for such a 
mass gap, Sm 2 ~ 3 x 10~ 3 eV 2 , in view of the large mixing angle, the sterile neutrino would 
come into equilibrium with the particle thermal bath via — v s oscillations 0. 

In this paper we formulate a mechanism capable of blocking the active-sterile neutrino 
oscillations that operates in the framework of theories where the lepton number is spon- 
taneously broken. It has been shown |18j that a generic feature of these theories is the 
production of coherent, long-range Majoron fields. While they can appear in stars as a 
result of neutrino oscillations or any other lepton number violating large scale process, we 
found that in the Early Universe a cosmological Majoron field emerges also as a result of 
a primordial lepton asymmetry carried by the scalar particles and complex scalar field (a) 
whose expectation value spontaneously breaks the lepton number. The Majoron field am- 
plitude is thus naturally proportional to the lepton number of the Universe, and baryon 
number as well, due to the L and B violating sphaleron processes. 

The leptons, and neutrinos in particular, interact with the derivatives of the Majoron 
field, which gives rise to new neutrino potentials that are relevant for the oscillation phe- 
nomena. The potentials are inversely proportional to the second power of the scale of lepton 
number symmetry breaking ((a)) but this scale can be much smaller than the electroweak 
breaking scale. For a lepton number breaking scale of the order of 1 GeV, a Majoron field 
associated with a lepton number asymmetry of the same order of magnitude as the baryon 
asymmetry can block the active-sterile neutrino oscillations at temperatures in the MeV 
range for a neutrino mass gap Sm 2 ~ 3 x 10~ 3 eV 2 . 

This provides an interesting way to block the sterile neutrino oscillations as it does not 
require the extraordinarly high lepton asymmetries, ~ 1CU 5 , that are necessary || in the 
framework of the standard model weak interactions. In fact, the standard model neutrino 
potentials that are proportional to the background lepton asymmetries are suppressed by 
the Fermi constant. In contrast, the potentials due to a Majoron field vary as the inverse 
square lepton number symmetry breaking scale. For a 1 GeV energy scale one immediately 
obtains the five orders of magnitude increase factor that brings 10 -5 down to the baryon 
number asymmetry B ~ 10 -10 . 
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